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We propose a stochastic model for the intersection of two urban streets. The vehicular traffic 
at the intersection is controlled by a set of traffic lights which can be operated subject to fix-time 
as well as traffic adaptive schemes. Vehicular dynamics is simulated within the framework of the 
probabilistic cellular automata and the delay experienced by the traffic at each individual street is 
evaluated for specified time intervals. Minimising the total delay of both streets gives rise to the 
optimum signalisation of traffic lights. We propose some traffic responsive signalisation algorithms 
which are based on the concept of cut-off queue length and cut-off density. 
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I. INTRODUCTION 



Modelled as a system of interacting particles driven 
far from equilibrium, vehicular trajfic provides the pos- 
sibility to study various aspects of truly non-equilibrium 
systems which are of current interest in statistical 
physics^^^. For almost half century, physicists have been 
challenged to understand the fundamental principles gov- 
erning the vehicular flow'*^''. Recently, discrete mod- 
els such as cellular automata have provided a signif- 
icant theoretical framework for the discipline of traffic 
flow modelling. The first cellular automata, the so-called 
BML model, was introduced by Biham, Middleton and 
Levine. This model described a simplified network of ur- 
ban intersections®, soon after, cellular automata, found 
their way in highway traffic through the pioneering work 
of Nagel and Schreckenberif which became the ancestor 
of many papers in the literature ( for a review see ref.^'^ 
and the references therein ). The BML model itself was 
later generalized to take into account several realistic fea- 
tures such as faulty traffic lights, independent turning 
of vehicles, and green- wave synchronizations'^"^''. The 
Nagel- Schreckenberg and BML model were recently com- 
bined to cast an upgraded version of city traffic models^®. 
In a very recent paper, the model is now extended to ac- 
count for different types of global signalisation^^. 
Despite the aforesaid efforts and those carried out by traf- 
fic engineers, the subject of optimal signalisation of re- 
alistic urban networks has not yet been comprehensively 
reviewed. In the above approaches, the main concern has 
been focused on the global strategies of the traffic net- 
work and frequently the role of isolated intersection have 
been suppressed. We believe that the optimisation of 
traffic flow at a single intersection is a substantial ingredi- 
ent towards a global optimisation. Isolated intersections 
are fundamental operating units of the sophisticated and 
correlated urban network and thorough analysis of them 
would be advantageous toward the ultimate task of the 
global optimisation of the city network. In this regard, 
our objective in this paper is to analyze the traffic state 



of an isolated intersection in order to find a better insight 
into the problem. In addition to theoretical viewpoint, an 
investigation of isolated intersection could be of practical 
importance. To a very good approximation, marginal in- 
tersections in cities are unaffected by other intersections 
and can be regarded as isolated ones. Generally there 
are two basic types of control for traffic lights at inter- 
sections: fixed-cycle and traffic-responsive^'^'^^ . Both of 
these methods can be implemented via centralized or de- 
centralized strategies. The application of each method 
strongly depends on traffic condition and the topology of 
the city network. In this paper we study the basic fea- 
tures of traffic flow at a single intersection which is con- 
trolled both under fixed-cycle as well as traffic responsive 
schemes. 



II. FORMULATION OF THE MODEL 

An isolated intersection is formed at the intersection of 
two streets. The streets, in principle, can each carry two 
opposite flows of vehicles. Depending on the designing 
of the intersection, different phases of movement can be 
defined ( a phase of traffic is defined as the fiow of vehi- 
cles that proceed an intersection without conflict). Here 
for simplicity we restrict ourselves to the simplest struc- 
ture: a one-way to one-way intersection. With no loss 
of generality, we take the direction of the flow in the first 
street, hereafter referred to as the street A, northwards. 
The other street ( hereafter referred to as street B) con- 
ducts a one-way eastward fiow. Cars arrive at the south 
and west entrances of the intersection. Space and time 
are discretized in such a way that each chain is divided 
into cells which are the same size as a typical car length. 
Time is assumed to elapse in discrete steps. We take the 
number of cells to be L for both roads. Each cell can be 
either occupied by a car or being empty. Moreover, each 
car can take discrete- valued velocities 1, 2, • ■ • , Vmax- To 
be more specific, at each step of time, the system is char- 
acterized by the position and velocity configurations of 
cars and the traffic light state at each road. 
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The system evolves under a generalized discrete- 
time Nagel-Schreckenberg (NS) dynamics. The gener- 
alized model incorporates anticipation effects of driving 
habits^^. ft modifies the standard NS model at its second 
step i.e., adjusting the velocities according to the space 
gap. Let us briefly explain the updating rules which are 
synchronously applied to all the vehicles. We denote the 
position, velocity and space gap (distance to its leading 
car) of a typical car at discrete time t by x^^^ , ti*-*-* and g*^*^ 
and the same quantities for its leading car by xf' , vf* and 

gf' . Assuming that the expected velocity of the leading 
car, anticipated by its follower, in the next time step t + l 
takes the form i'[*^„fj = min{g\^\v\^'')^ we define the ef- 
fective gap as g'^^^jf := g^'^^ +max{vl^l^^.- gapsecure, 0) in 
which gapsecure is the minimal security gap. Concerning 
the above considerations, the following updating steps 
evolves the position and the velocity of each car. 

1) Acceleration: 

2) Velocity adjustment : 

3) Random breaking with probability p: 

if random < p then := maxiv'^*'^'^/'^^ — 1, 0) 

4) Movement : := a;^*) + 

Let us now specify the physical value of our time and 
space units. Ignoring the possibility of existence of long 
vehicles such as buses, trucks etc, the length of each 
cell is taken to be as 5.6 metres which is the typical 
bumper-to-bumper distance of cars in a waiting queue. 
Concerning the fact that in most of urban areas a speed- 
limit of 60 kilometre/hour should be kept by the drivers. 



we quantify the time step in such a way that Vmax = 6 
corresponds to the speed-limit value (taken as 60 km/h). 
In this regard, each time step equals two seconds and 
therefore each discrete increments of velocity signifies a 
value of lOkm/h which is equivalent to a comfortable 
acceleration of 1.4 rajs^. We have also set the horizon 
length L = 70 cells and gap sec = 1- The state of the 
system at time < -I- 1 is updated from that in time t by 
applying the modified NS dynamical rules. 

Step 1 : signal determination. 

We first specify the signal states for all of the driving 
directions. In subsequent sections we will, in detail, ex- 
plain the scheme at which the traffic lights change their 
colour. 

Step 2 : movement in the green road. 

At this stage, we update the position and velocities of 
cars on the green road according to the movement rules 
which are synchronously applied to each car . 

Step 3 : movement at the red road, delay eval- 
uation . 

Here the updating is divided into two parts. In the 
first part, we evaluate the delay of cars waiting on the 
red period of the signalisation. In the second half, we 
update the position and velocities of the moving cars 
approaching the waiting queue. We should note that 
once the signal switches to red, the moving cars continue 
their movement until they come to a complete stop when 
reaching to the end of the waiting queue. As soon as a 
car comes to halt, it contributes to the total delay. In 
order to evaluate the delay, we measure the queue length 
(the number of stopped cars) at time step t and denote 
it by the variable Q. We recall that pos'''^'^[i, t] — 1 for 
i — 1- ■ - Q and zero at i = Q + 1. Delay at time step < + 1 
is obtained by adding the queue length Q to the delay at 
time step t. 

delay{t+l) = delay{t) + Q{t) (1) 

This ensures that during the next time step, all the 
stopped cars contribute one step of time to the delay. 
The next part describes the update of positions and ve- 
locities of moving cars. Moving cars can potentially be 
found in the cells Q + 2,Q + 3, ■ ■ ■ , L. We update their 
positions and velocities accordingly. 

step 4 : entrance of cars to the intersection. 

So far, we have dealt with those cars within the hori- 
zon of the intersection which goes up to the entry point 
located at site L. Here we discuss the entrance of cars 
into the intersection. Generally speaking, during each 
green period, a fraction of the queue will dissolve and 
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pass the intersection. If the average arrival rate of cars 
exceeds the maximum evacuation rate of the lane, then, 
on an average level, a fraction of a generic queue will 
not be able to go through the intersection and should 
wait until the next upcoming green period arrives. Cor- 
respondingly in the course of time, the remainders accu- 
mulate giving rise to a growing queue length. In reality 
we rarely observe this phenomena since an actual inter- 
section is linked to other ones hence the possibility of such 
a rare event is restricted to very exceptional cases where 
there is an overwhelming large amount of incoming flux. 
Throughout the paper, we assume that the average in- 
flow rate is sufficiently below this evacuation rate so that 
a generic queue will have enough green time to dissolve 
and that the intersection is able to support the incoming 
flux. In this case, there is a typical queue length which 
affects the motion of incoming cars. The motion of a car 
approaching the red light is affected via two factors: the 
distance to the end of the queue and the traffic light sig- 
nal which act as a hindrance to the incoming traffic flow. 
In the light traffic state under consideration, there exist 
an interaction distance to the queue beyond which one 
can, to a good extent, assume that the cars are moving 
without being affected by the traffic light signal. Away 
from this interaction zone, the cars move according to 
the movement rules without any hindrance. We take the 
position of the place at which the cars enters the horizon 
of the intersection to be 70 cell equivalent to 400 metres. 
The time head-ways between entering cars at this entry 
location vary in a random manner which consequently 
implies a random distance headway between successive 
entering cars. As a candidate for describing the statisti- 
cal behaviour of random space gap of entering cars, we 
have chosen Poisson distribution. The Poisson distribu- 
tion function has been used in a variety of phenomena 
incorporating the modelling of " queue theories" and has 
proven to be a good estimation of reality^'^. In addition, 
it has the merit of taking only discrete values which is de- 
sirable to us in the view of the fact that in our model the 
gap is a discrete variable. According to this distribution 
function the probability that the space gap between the 
car entering the intersection horizon and its predecessor 
be n is : p{n) = — where the parameter A speci- 
fles the average as well as the variance of distribution 
function. The parameter A is a direct measurement of 
traffic volume. A large value of A describes a light traffic 
while on the other hand, a small-valued A corresponds 
to a heavy traffic state. Now wc let the cars enter into 
the intersection's horizon. At the end of the movement 
rule, we evaluate the position of farthest cars on both 
streets. We denote them by last a and lasts- By defini- 
tion , posA[j,t] = for j > last A. A similar statement 
applies to street B. In order to simulate the entrance of 
cars into the horizon of the intersection, we randomly 
choose an integer weighted by Poisson distribution func- 
tion. This number represents the gap of the oncoming 
successor of the farthest car. Let us denote these num- 



bers (headways) by Ha and hs for street A and B re- 
spectively. Once the random gap is chosen, we create a 
car at the position lastA + Iia {lasts + hs) of the street 
A (B) respectively. The created car survives provided 
the following constraint is satisfied: last a + Ha < L 
( lasts + hs < L). If the position of the created car 
exceeds the horizon length L, the creation procedure is 
rejected. The above ad hoc rules updates the configura- 
tion of the intersection in time. In the next section, we 
will explain our simulation results. 



III. SIGNALISATION OF TRAFFIC LIGHTS: 
FIXED TIME SCHEME 

In this controlling scheme, the traffic flow is controlled 
by a set of traffic lights which are operated in a fixed- 
cycle manner. The lights periodically go green with a 
fixed period (cycle length) T. This period is divided into 
two parts : in the first part, the traffic light is green for 
street A (simultaneously red for street B). This part 
lasts for Tg seconds { Tg < T). In the second part, the 
lights change colour and the movement is allowed for 
the cars of road B. The second part lasts from Tg to 
T. This behaviour is repeated periodically. Cars enter 
the intersection and a fraction of them experience the 
red light and consequently have to wait until they are 
allowed to go through intersection during the upcoming 
green periods. Now we raise the basic question: 

how should one adjust the ratio of s^sil {ji order to 
optimise the throughput flow?. 

There is now an almost well-established agreement on 

the quantitative definition of optimisation. Borrowing 
from the traffic engineering literature, we adopt opti- 
mised traffic as a state in which the total delay of ve- 
hicles is minimum. In one of our preceding works^'', we 
analytically evaluated the total delay in terms of arrival 
as well as that of the exit rates of vehicles. However, 
our approached was based on the simple assumption of 
the time-constancy of the arrival rates. In reality, we 
know that successive cars arrive with fluctuating time- 
headways which consequently induces time-varying ar- 
rival rates. In this paper we address the question of non- 
constant rates. In order to evaluate the delay, we have 
simulated the flow of vehicles. For the sake of simplic- 
ity, we have assumed that each street has a single lane. 
For streets with more than one lane, one simply should 
multiply the value of delay by the number of lanes. 



A. Simulation Results: Symmetric inflow 

Wc let the intersection evolve for 1800 time steps which 
is equal to a real time period of one hour. We let the 
green time of street A, Tg, vary from zero to T. For 
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each value of Tg , we evaluate the aggregate delay corre- 
sponding to both traffic lights during 1800 time steps. 
We have averaged the results of 50 independent run of 
the programme. Let us first consider the symmetric traf- 
fic states in which the traffic conditions are equal for 
both roads. In this case, we equally load the intersection 
with entering cars spatially separated by random space 
gap (Poisson statistics) from each other. The following 
graphs depicts the total delay curves as a function of Tg 
allocated to road A for two values of cycle length. 
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Fig. 2 and 3: The cycle length is 40 seconds (top graph) 
and 60 seconds (bottom graph). Total delay versus the green 
time of road A is sketched for various in-flow rates. 



The general shape of the total delay curve resembles 
a "U-shaped" valley. The middle area, where the total 
waiting time is minimised, corresponds to a situation in 
which the evacuation rate of roads exceeds the in-flow 
and the queues can be dissolved during one green pe- 



riod and consequently all the waiting cars will be able 
go through the intersection in the upcoming green pe- 
riod. The optimal traffic flow is obtained by keeping the 
lights at equal timing Tg = -j. Before preceding further, 
it would be useful to discuss, in detail, the conditions in 
which the intersection would not able to support the in- 
flux and hence the queue length starts growing. Suppos- 
ing the fraction ^ of the cycle is given as green time to a 
street say A. Therefore during the one hour time interval, 
street A receives 3600^ seconds of green times. Using 
the fact that in green phase vehicles are going through 
the crossing with the approximate rate of 0.5/sec/lane, 
on an average level, the maximum out-flow capacity, de- 
noted by (Cmax) would be: {Cmax) = 1800^. On the 
other hand, total number of in-flow can be estimated 
once the average space gap is given. To do so, we first 
obtain the average time headway of entering vehicle (in 
second), denoted by Aa_4, as follows Aa_4 — 5.6 ^'^'^ . 
In the above formula, we have assumed that entering 
cars have the maximum velocity and 5.6 denotes the cell 
length in metre. Concerning the above consideration, 
during one hour, the average number of entering cars are 
approximated by {Nin) — -P^o . Avoiding the occurrence 
of growing queue, leads to the satisfaction of the following 
constraints^: 



{Nin) !i {Cmax) 



(2) 



Substituting maximum velocity by 16.7 m/s in the 
above constraint, gives rise to condition Tg > 1^. Taking 
T = 60 seconds, the minimum consistent value of Tg for 
A — 13, 16 and 19 cells would be 27, 22 and 19 seconds 
respectively. Alternatively, the lowest \a is restricted 
to 12 cells or 67 metres. Similar arguments should be 
applied to street B. Concerning the fact that the green 
time of street B is T — Tg, one simply deduce that the 
condition T — Tg > should hold. Combining the two 
conditions on Tg, one arrives at the following inequality 
on Tg for a stationary condition of the queues. 



6T . 6 , 

T-<Tg<Til- — ) 



(3) 



Consequently, the allowed Xa, and Tg should sat- 
isfy the above inequality, otherwise the queue will grow 
infinitely in time. From the above relation one con- 
cludes two necessary conditions on the rates: As > 6 
and t'- -I- -r- < ^ which hold regardless of the value of 
green time. We note that the above arguments are based 
on simple mean-field approximation and no fluctuation 
has been taken into account. There are two origins of 
fluctuation, the first one concerns the stochastic space 
gap of cars which make the incoming flux deviated from 
a constant one, and second, the interaction among cars 
and randomness arising from the car movements. There- 
fore the region close to critical volume A^ = 12 needs 
special treatment. Simulation results shows that away 
from the critical region, fluctuations do not violate our 
mean-field conclusions. 
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B. Asymmetric in- flow rates 



In the asymmetric states, the entering rates into the 
roads differ from each other. An interesting asymmetric 
state is the intersection of a major to a minor street. A 
large fraction of urban intersections Ues in the category 
of major-to-minor. The signahsation of these types of 
crossing is still a controversial subject. The main reason 
is that in most of intelligent real-time controller systems, 
the signalisation of these intersections are highly affected 
by the major intersections signalisation schemes which 
frequently overlook the local optimisation of minor in- 
tersections. In our model, a major-to-minor crossing is 
modelled by a light traffic in one road and a congested 
one in the other road. The following graph depicts the 
behaviour of the delay curve at a major-to-minor inter- 
section. As observed, the minima of delay curves are 
shifted toward right which expresses the fact that the 
majority of green times should be given to the road with 
higher in-flow rate. Analytical considerations lead to 
distribution of green times in proportion to the in-flow 
rates of the roads^"'. To be more specific, let us denote 
the average arrival rate at street A and B by aA and 
as- Note that a a and as are proportional to inverse of 
entering space gaps. We now evaluate the waiting time 
of street A during one complete cycle. The number of 
cars arriving the queue during the time interval [t, t -f dt] 
in the red period, which lasts T ~ Tg seconds, is aAdt. 
These cars should wait T — Tg — t seconds until the up- 
coming green period starts. Therefore, their contribution 
to delay reads aAdt{T — Tg — t). The aggregate delay of 
street A during the cycle is obtained by summation of 
infinitesimal delays as follows: 

aAdt{T -Tg~t)^ iaA(T - Tgf (4) 

Similarly the contribution of street B to total delay can 
simply be evaluated as ^asTg. Adding the street delays 
together and minimising it with respect to Tg gives rise 
to the optimal Tg-. 



rj^opt ^ '±A rj. (-5^ 

aA + as 

The above result corresponds to distribution of green 
times in the ratio of in-fiow rates. Our simulations re- 
sults are in good agreement to this conclusion. In fact 
the position of optimal green times, i.e., the minima of 
the delay curves, grow linearly with the linear decrease of 
in-fiow rates which are in turn proportional to the inverse 
of average space gap. 
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Fig. 4: Total delay for asymmetric traffic volumes in a 
major-to-minor crossroads. The space gap of the major road 
is set to 13. The space gap of the minor road is varied. T=60 
seconds. 



IV. SIGNALISATION OF TRAFFIC LIGHTS: 
TRAFFIC ADAPTIVE 

We now discuss the traffic adaptive controlling scheme 
in which the light signalisation is adapted to the traffic 
at the intersection. Nowadays advanced traffic control 
systems anticipate the traffic approaching intersections. 
Traffic-responsive methods have shown a very good per- 
formance in controlling the city traffic and now a variety 
of schemes exists in the literature^"'^^'^^'^^. In these 
schemes, the data obtain via traffic detectors installed 
at the intersection is gathered for each movement di- 
rection and it is possible to count the queue lengths 
formed behind the red lights. One can also measure 
the time-headways between successive cars passing each 
lane detector. Thus it is possible to estimate the traffic 
volume existing at the intersection. There are various 
methods for distribution of green times. In what follows 
we try to explain some standard ones. In each scheme, 
the green time of a typical green street is terminated if 
some conditions are fulfilled. By green (red) street we 
mean the street for which the traffic light is green (red). 
We now state the termination conditions for each scheme. 

Scheme (1): The queue length in the confiicting di- 
rection exceeds a cut-off value Lc- This scheme only 
concerns the traffic states in the red street. 

Scheme (2): The global car density in the green street 
falls below the cut-off value pc- Here the algorithm only 
considers the traffic state in the green street. 
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Scheme (3) : Each direction is endowed with two con- 
trol parameters Lc and pc. The green phase is terminated 
if the conditions: < pc and L'' > Lc are both satisfied. 

Here the algorithm implements the traffic states in 
both streets. The superscripts "r" and "g" refer to words 
"red" and "green" respectively. We note that the first 
two schemes are special cases of the more general scheme 
(3). Schemes (1) and (2) are the limiting cases of schemes 
(3) by letting pc ^ I and Lc — > respectively. In gen- 
eral, the numerical value of control parameters and 
Pc could be taken different for each individual street. In 
what follows we present our simulation results for differ- 
ent types of signalisation schemes introduces above. 

A. Simulation results: symmetric inflow 

We let the intersection evolve for 1800 time steps which 
is equal to a real time period of one hour. We evaluate 
the aggregate delay for both streets as well the number 
of passed vehicle during 1800 time steps at each direction 
for different traffic situations. As the first situation, we 
have considered the symmetric traffic state in which the 
traffic conditions are equal for two streets. In this case, 
we equally load the streets with entering cars spatially 
separated by random space gap, obeying Poisson statis- 
tics, from each other. We first discuss the first scheme 
in which the lights change colour as soon as the queue 
length reaches the cut-ofi' length. We set an equal cut- 
off length for both streets. The following graph depicts 
the total delay curve as a function of traffic volume for 
various cut-off lengths. 



Fig. 5: total delay in terms of average space gap of entering 
cars for some cut-off queue length Lc — 5, 6, 7, 8. 

We observe that for each cut-off length, total de- 
lay shows a slight decrease with respect to the in-ffow 
rate. The graph demonstrate that the shorter the cut-off 



length, the less the total delay. The next graph shows the 
delay versus in-flow rate when the lights are signalised ac- 
cording to scheme (2). Here the green light is terminated 
below a certain occupancy in the corresponding moving 
ffow. 
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Fig. 6: total delay in terms of average space gap of entering 
cars for some cut-off densities. 

We observe that the first scheme gives rise to lower 
delay. The reason is that in relatively low traffic volume, 
more cars have to stop in the red light in order that the 
global density in the green direction reaches to its cut-off 
value. This raises the delay compared to scheme (1). 
Finally we discuss the third scheme in which two condi- 
tions should be satisfied for termination of green light. 
In this scheme both cut-off densities and cut-off lengths 
are taking into account. We have investigated the case 
where cut-off parameters are equal for both streets i.e., 
Pc=Pc= Pc and l[^'^ = Lf ^ = Lc. The following 
table shows the delay table of one-hour performance for 
various values of pc and Lc- Traffic in-flow rates are 
= As = 13. 



Pc i \ Lc^ 


5 


6 


7 


8 


9 


0.10 


14000 


15500 


15900 


17900 


18600 


0.20 


13200 


14900 


14700 


16500 


17400 


0.30 


11500 


12800 


13200 


15500 


16400 


0.40 


11200 


11900 


12600 


14000 


15300 


0.50 


10900 


11200 


11700 


12700 


13900 



The above result sets the optimal parameters at pc — 
0.5 and ic = 5. We have tested a variety of other traf- 
fic volumes. The results demonstrate that in the control 
parameter space, the optimal region lies in high pc and 
low Lc- The above result shows that scheme (3) gives an 
improved delay with respect to scheme (2). Our simula- 
tion results shows that this conclusion can be made for 
other values of A hence the third scheme is more opti- 
mal than the second one. However the results show that 
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still the first scheme gives lower delays compared scheme 
(3). Specifically for A = 20 cells, the predicted delays 
are 8500, 10100 and 9600 seconds respectively. The next 
graphs depict the behaviour of consecutive green times 
of streets. Wc have chosen the green times of street A. 
Fluctuating traffic volume at streets induces fluctuating 
green periods which are known as the green time plans 
in the traffic engineering literature^^'^^. We note that a 
realistic optimising algorithm must satisfy the constraint 
that the duration of green phases should not be shorter 
that a minimal amount e.g. 10 sec or so. This mini- 
mal value reflects the actual time required for an immo- 
bile waiting queue to accelerate and make a considerable 
movement forward. Our simulation results shows that 
Lc > 4 corresponds to typical green periods greater than 
12 seconds. It may be useful to analyze the statistics 
of green plans. Denoting the consecutive green times 
by T^'^\t^'^\t^^^ • • •, We first obtain the basic statisti- 
cal properties which are the moments of the distribution 
of the green times. The bottom figure shows the green 
plan histogram of street A for = = 13. For the 
mentioned traffic state, the average green time and the 
standard deviation of street A are 18 and 2.54 seconds 
respectively. 
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Fig. 7: Green times liistogram of street A. Aa = As = 13 
cells. The lights are controlled via scheme (1). The cut-ofT 
length is taken as 7 cars. 



B. Asymmetric in- flow 

scheme (1): 

We now investigate the situations where the streets are 
not equally loaded. For this purpose, we fix the in-flow 
rate of street A at Aa = 13 and vary the in-flow rate of 
street B. The following graph exhibits the total waiting 
time curves as well as those of each street for some cut- 



off lengths (taken equal for both streets) . Here again the 
less delay is achieved for shorter cut-off lengths which in 
turn correspond to shorter green phases. 



Fig. 8: Total delay versus the inverse traffic volume of 
street B for various cut-off lengths. \a is fixed at 13. 

We observe that decreasing the traffic volume in street 
B leads to a decrease in total delay as expected. In the 
following figure, the dependence of each street delay on 
the traffic volume of street B is exhibited in details. One 
observes that the delay of street B starts growing whereas 
that of street A decreases. This shows that the algorithm 
does not always act optimally for individual directions. 
In fact, the delay of street A, which is more congested 
than street B, behaves decreasingly while the delay of 
the less congested street, B, behaves increasingly. Nev- 
ertheless, the algorithm acts in an optimal manner when 
taking into account the whole intersection. 



Fig. 9: Delay of individual streets as a function of inverse 
traffic volume of street B. Traffic volume of street A is kept 
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fix Xa ~ 13. 



for various controlling densities. 



We have also examined the cases where cut-off lengths 
take different values for streets with non-equal incom- 
ing flux. Simulation results shows that optimal cut-off 
lengths should be taken equal and as short as possible. 
For instance, in the case Xa = 14 and As = 24 the min- 
imised total delay when Li"^' = = 5 is 8500 seconds. 

scheme (2); 

Let us now investigate the predictions of scheme (2). 
As explained, in this method, once the global density of 
the moving direction (in the vicinity of the intersection) 
falls below the cut-off density, green phase is terminated 
irrespective of the queue length in the red direction. The 
global density is obtained by measuring the occupancy of 
Lp cells before the crossing point. The following graphs 
shows the one-hour delay behaviour of individual streets 
as well as that of the entire intersection for various con- 
trolling densities. Analogous to the predictions of first 
scheme, here the algorithm increases the delay in the mi- 
nor road whereas the delay is decreased in the major road 
upon decreasing the traffic volume in the minor road. 



cut-off density= 0.5 , street A 
cut-off density= 0.5 , street B 
cut-off density= 0.4 , street A 
cut-off density= 0.4 , street B 




Average space gap of street B (cells) 



cut-off density= 0.2 
cut-off density= 0.3 
cut-off density= 0.4 
cut-off density= 0.5 




Average space gap of street B (cells) 



One observes that the delay predictions of the second 
scheme are higher than those of the first scheme. We have 
also examined the cases where cut-off densities of each 
street take non-equal values. The following table exhibits 
the values of one-hour total delay for various amounts 
of cut-off densities of each street. Traffic volumes cor- 
respond to = 13 and As = 18 and Lp is set to 10 cells. 
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0.30 


12900 


13300 
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0.40 


11400 


12100 
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14200 


14000 


0.50 


10800 
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11800 


12400 
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0.60 


10600 


10800 


10700 


10900 


11000 



We observe that optimal flow corresponds to set 
PA — 0.1 and pb = 0.6. In contrast to scheme (1) 
where optimal cut-off lengths are equal to each other, 
here we see that in asymmetric traffic volume, the opti- 
mal cut-off densities are non-equal. We note that other 
values of cut-off densities give rise to unrealistic green 
time plans and hence are excluded from the table. It 
would be useful to compare the minimal delay with the 
predicted one in the scheme (1). For the traffic volumes 
corresponding to A^i = 13 and As — 18, scheme (2) gives 
the minimal delay at 10600 seconds whereas the scheme 
(1) predicts the value 8700 seconds for Lc = 5. We have 
also obtained the total delay table for another asymmet- 
ric traffic volume where in-flow rates are A^i = 13 and 
As = 25. For this choice of traffic volumes, the optimal 
cut-off densities are pA — 0.10 and ps — 0.40 which 
are relatively close to the preceding optimal values. The 
values of minimal delays are 9300 and 8500 seconds in 
scheme (2) and (1) respectively. We have investigated a 
variety of traffic states corresponding to different in-flow 
rates. The results show that scheme (1) gives better op- 
timal states than scheme (2). Before discussing scheme 
(3), we would like to discuss the role of parameter Lp. 
In fact, the measurement of global density depends on 
the length Lp. This length itself could be regarded as 
a control parameter. We have carried out simulations 
for various values of Lp. Our results shows that the 
Lp = 10 corresponds to the best choice. For instance the 
result obtained for the same traffic state investigated in 
the above table and Lp = 12 gives the optimal delay as 
11000 seconds which is higher than the value predicted 
by Lp = 10. 

scheme (3): 



Fig. 10: Top: delay of individual streets as a function of 
inverse traffic volume of the minor street B. Traffic volume of 
street A (major street) is kept fix Aa = 13 cells. Bottom: to- 
tal delay versus the inverse traffic volume in the minor street 



Finally we discuss the third scheme. In this scheme 
both cut-off densities and cut-off lengths are taking into 
account. The following graphs shows the one-hour de- 
lay behaviour of individual streets as well as that of the 
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entire intersection for various controlling densities. Anal- 
ogous to the predictions of first and the second schemes, 
here the algorithm increases the delay in the minor road 
whereas the delay is decreased in the major road upon 
decreasing the traffic volume in the minor road. We have 
taken the cut-off parameters equal for both streets i.e., 
= Lf and = pf . 



cut-off lengtfl= 6, street A 
cut-off lengtfl= 6, street B 
cut-off lengtfl= 8, street A 
cut-off lengtfi= 8, street B 




15 20 25 

Average space gap of street B (cells) 



cut-off lengths 5 cars 
cut-off lengths 6 cars 
cut-off lengths 7 cars 
cut-off lengths 8 cars 




15 20 25 

Average space gap of street B (cells) 



Fig. 11: Top: delay of individual streets as a function of 
traffic volume of the minor street B. Traffic volume of street 
A (major street) is kept fix Aa ~ 13 cells. Bottom: total 
delay versus the traffic volume in the minor street for various 
controlling densities. Cut-off densities are taken 0.5 for both 
streets. 

Another investigation considers the effect of cut-off 
densities for fixed values of cut-off lengths. The follow- 
ing table exhibits the delay when = 13 and Xb = 18. 
Lc = 5 while pc is varied. 
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The table above shows that optimal traffic ffow is ob- 
tained in the high-high region of iPc,Pc) sub-space of 
control parameters. The minimised delay is 9900 sec- 
onds. For the same in-flow rates, the values of minimal 
delays predicted by the schemes (2) and (1) are 10600 
and 8700 seconds respectively. This comparison shows 
that the third scheme acts better than scheme (2) but 
is less efficient than the first scheme. Further simulation 
results carried out for other traffic states demonstrates 
that for A_B < 19 (A^ fixed) scheme (3) acts better than 
scheme (2) but is less efficient than scheme (2) for higher 
in-flow rate of the minor street. Nevertheless, scheme (1) 
gives the least delay for all asymmetric states. 



V. SUMMARY AND CONCLUDING REMARKS 



In this paper we have investigated the optimisation of 
vehicular traffic ffow at an isolated signalised intersection 
in the framework of cellular automata. Marginal urban 
areas are places where single intersections are frequently 
designed and operated. Investigations on single junctions 
can be of practical relevance for various applications in 
city traffic. Basically there are two methods of signalisa- 
tion of the traffic lights: fixed-time or adaptive. In traffic 
adaptive strategies, it has always been a subject of argu- 
ment whether to control an intersection under a central- 
ized or decentralized scheme. In special circumstances, 
decentralized local adaptive strategies operate more effec- 
tively than globally adaptive strategies^°'^^'^^~^^ and of- 
ten show a very good performance. For this purpose, we 
have developed and analysed prescriptions for the traffic 
light signalisation at single intersections. Our simulation 
results confirm that in fixed-time scheme, green times 
should be distributed proportional to the traffic in-flow 
rates. We have simulated and analysed three traffic re- 
sponsive signalisation schemes. The signalisation mecha- 
nisms are based on the concepts of cut-off queue lengths 
and densities. Two major conclusion can be made from 
our simulation results. First, the traffic adaptive scheme 
act more optimally than the fixed-time. Secondly, the 
best traffic adaptive scheme is the one in which the flow 
in the green direction is terminated as soon as the queue 
length in the opposite direction exceeds the cut-off value 
Lc- In this method, the traffic states in the green di- 
rection are not taken into account. It should be noted 
that the results in this paper have been obtained un- 
der a modified version of Nagel-Schreckenberg model. 
Whether these result are robust against more realistic 
vehicle movement models needs more exploration and is 
the subject of our current investigation. Throughout the 
paper we have used the Poisson statistics for the head- 
way of approaching cars. However the realistic in-flow 
traffic would certainly deviate from this statistics. Inves- 
tigations on the impact of the entrance statistics on the 
above results will shed more lights onto the problem. 
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